Bubble dynamics in N dimensions 
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Cavitation and bubble dynamics are central concepts in engineering, the natural sciences, and 
the mathematics of fluid mechanics. Due to the nonlinear nature of their dynamics, the governing 
equations are not fully solvable. Here, the dynamics of a spherical bubble in an iV-dimensional fluid 
are discussed in the hope that examining bubble behavior in N dimensions will add insight to their 
behavior in three dimensions. Several canonical results in bubble dynamics are re-derived, including 
the Rayleigh collapse time, the Rayleigh-Plesset equation, and the Minnaert frequency. Numerical 
simulations are used to examine the onset of nonlinear behavior. Overall, the dynamics of bubbles 
are faster at higher dimensions, with nonlinear behavior occurring at lower amplitudes. Several 
features are found to be unique to three dimensions, including the trend of nonlinear behaviour and 
apparent coincidences in timescales. 



INTRODUCTION 

The behavior of bubbles and cavities has been an area 
of interest in fluid dynamics since the time of Leonardo 
da Vinci Serious theoretical research began in the 
early twentieth century when Lord Rayleigh derived their 
collapse time in order to study the effect of cavitation 
damage to ship propellers Q . Bubbles and cavitation re- 
main an active area of research in fluid dynamics [3| and 
soft condensed matter physics Notably, bubbles find 
application as contrast agents for biomedical ultrasound, 
where their oscillations can be detected acoustically [5j. 
The mechanism of sonoluminescence, a flash of light re- 
leased by a collapsing bubble, is still an unsolved problem 
in physics Q. 

The dynamics of gas-filled spherical bubbles are de- 
scribed using the Rayleigh-Plesset equation, a special 
case of the Navier-Stokes equation that describes the os- 
cillation of a spherical cavity in an infinite incompressible 
fluid. In the linear regime the bubble oscillates harmoni- 
cally, pulsating with a breathing mode given by the char- 
acteristic Minnaert frequency. At higher oscillation am- 
plitudes, the bubble exhibits strong nonlinear behavior. 
Like the Navier-Stokes equation, there is no known gen- 
eral solution to the Rayleigh-Plesset equation or to the 
Rayleigh-like collapse of a bubble. The equations are of- 
ten solved numerically, but in many cases the oscillations 
may become chaotic. 

Higher or lower dimensional analyses are often utilized 
to aid understanding of three dimensional physical prob- 
lems. For example, the Ising model of ferromagnetism, 
while not solvable in three dimensions, can be exactly 
solved in two dimensions and is solvable with mean-field 
theory in more than four dimensions 0] ■ Varying the di- 
mensionality of the system as a free parameter can often 
provide fresh physical or mathematical insight [8] . Higher 
dimensional generalizations have previously been utilized 



to solve theoretical problems in acoustics: the question 
"Can one hear the shape of a drum?" was first answered 
on a sixteen dimensional torus 9]. Regularity criteria 
established for the Navier-Stokes equations in N dimen- 
sions [1(3] have been used to establish such criteria for 
the three-dimensional case ll|. Bubble dynamics, par- 



ticularly nucleation and expansion, are also of interest in 
high-energy physics in the areas of relativistic hydrody- 
namics [121 ] and cosmological phase transitions in three 



or more spatial dimensions 13]. The stability of fluid 



droplets and cylinders were analyzed in N dimensions 
by Cardoso and Gaultieri 14[ and differential equations 
resembling certain types of iV-dimensional bubble were 
analyzed by Lima et al. [l5j, but their analysis focused 
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on the three-dimensional case. Notably, Prosperetti dis- 
cussed Rayleigh-like and oscillatory bubble dynamics in 
N dimensions but only focused on the three-dimensional 
case and presented inconsistent results [H. 

In this paper, the dynamics of a hyperspherical bubble 
are analyzed. Relevant terms such as the Rayleigh col- 
lapse time, Rayleigh-Plesset equation, and the Minnaert 
frequency are derived. A recent paper by Obreschkow 
et al. 0] presenting analytical approximations to bubble 
collapse will also be analyzed in its generalized form. As 
a case study into higher dimensional bubble dynamics, 
the onset of nonlinear behavior is analyzed. Examining 
the behavior of bubbles as a function of dimensionality, it 
is hoped that higher dimensional bubble dynamics cast 
light on unknown aspects of three-dimensional bubble 
dynamics. 

The derivations in this paper follow those historically 
used for three dimensions. Various assumptions are 
made: the bubble remains spherical at all times, the fluid 
is incompressible and infinite, the gas inside the bubble 
is an ideal gas with uniform pressure and constant num- 
ber. When comparing analytical results to numerical so- 
lutions, the Runge-Kutta method was used for second- 
order differential equations, and the Euler method was 
used for first-order differential equations. 
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DERIVATIONS 
The Continuity Equation in N Dimensions 

Before the dynamics of bubbles are re-derived in ar- 
bitrary dimension, it is necessary to understand the be- 
havior of fluids, and the associated vector calculus, in 
arbitrary dimension. Central to this is the continuity 
equation that formally states the conservation of mass, 
relating the velocity field u and the density p in a fluid: 



| + v-(H = o 



(i) 



In an incompressible fluid, p is constant and the equation 
is simplified: 



V • u = 



(2) 



In three-dimensional spherical coordinates with no angu- 
lar dependence, the divergence operator takes the form: 



t-, 1 d r 2 \ 

V • u = -it— (r u) 

r^ dr 



(3) 



Solving the continuity equation with spherical symmetry 
shows that the velocity field around a spherical source 
has an inverse-square fall-off: 



u{r) = -2 



(4) 



In N dimensions, the spherically symmetric divergence 
operator takes the form: 



1 d 



dr 



„n-i 



(5) 



and the spherical solution of the continuity equation is: 



u(r) = 



.N-l 



(0) 



This relationship will be used when deriving the kinetic 
energy of the fluid around a pulsating hyperspherical 
bubble. It is also worthwhile to express the surface area 
and volume of the unit sphere in N dimensions: 



SN = T(Nj2) Vn 



T N/2 



T(N/2 + 1) 



Rayleigh-like Collapse 



Rayleigh in 1917 Q derived a formula for the collapse 
time of a spherical cavity of radius R in an infinite liquid 
with pressure P and density p, considering the time it 
takes for the liquid to fill the cavity. The derivation lay 
in equating the kinetic energy of the fluid with the work 
performed by the bubble. The result was: 



^collapse — 0.91468 ■ R \j ~p 



(8) 



The collapse time can be derived in N dimensions: the 
kinetic energy is calculated by integrating along the ra- 
dial coordinate r the energies of each spherical volume 
element of density p and thickness dr as it expands or 
contracts in phase with the central bubble of radius R. 



KE = 



S N r 



N-l 



R 



R 



N-l 



dr = 



Snp R R 



2 D» 



2 N-2 
(9) 

The term in the square brackets represents the 
continuity-dependent velocity dispersion, ©. In two di- 
mensions, the kinetic energy integral does not converge; 
the dimensionality discussed in this paper is limited to 
N > 3. In the case described by Rayleigh, the pressure P 
in the fluid is constant during collapse. Thus, the work 
performed by the fluid is merely the change in volume 
multiplied by the pressure: 



W = V N P (R? - R N ) 



(10) 



Recalling that the ratio Vn/Sn is simply equal to 1/N, 
equating (TiTj|) and © and isolating R yields a first-order 
differential equation between R and t. 



dR 
~~dl 



P2(N 



P 



N 



nN 
Rp 

R N 



(11) 



The Rayleigh collapse time i_Rc, the amount of time re- 
quired for a bubble initially at R a to collapse to zero, can 
be found by isolating dt from (jTTJ) and integrating by dR 
from to R n : 



tRC 



p N 



R 



N 



P 2 (N - 2) J V R o -R N 



dR (12) 



Following the procedure used for TV = 3 by Rayleigh, the 
substitution /3 = is made: 



p N f 1 B N / 2 

to-^Tm^Lwy*" (13) 



N ■ Sn (7) This integral can be solved in terms of gamma functions: 




tRC — 



We define the unitless parameter r to simplify notation. 
For N = 3, r reduces to the familiar result of 0.91468, 
and is a monotone decreasing function as a function of 
N for N > 3. For large N, the gamma functions in 
the numerator and denominator simplify to ^/tt and 1 
respectively. The collapse behavior will be discussed in 
greater detail subsequently, but the chief result of this 
derivation is that the collapse time decreases in higher 
dimensions. 
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Analytical Approximations 

A recent paper by Obreschkow et al. Q examined 
Ray leigh- like collapse in greater detail. They expressed 
the collapse of a spherical bubble in normalized form, 
and derived a simple function to approximate it. Their 
approximation function, even at lowest-order, showed ex- 
cellent agreement with both the numerical solution of the 
Rayleigh equation and their experimental measurements. 
Here, the normalization and approximation will be gen- 
eralized for higher dimensions in order to fully explore 
the collapsing bubble behavior. 

The Rayleigh collapse equation is normalized with the 
variables r = R/R and t = t/tnc- It is expressed in a 
lesser known form relating the radius to its second deriva- 
tive: 

r = -T 2 r~ 4 (15) 

The normalized equation is generalized for arbitrary di- 
mension by differentiating (fTTT) with respect to t: 

r = -{N -2)T 2 r-^ N+ V (16) 

This allows an examination of the behavior of the col- 
lapsing radius for various dimensions independent of 
timescale. It can be seen that for higher dimensions, 
the bulk of the collapse occurs towards the end of the 
process, while at lower dimensions it is a more gradual 
process (Figure [1]). This is noted by Prosperetti [1|, who 
states that the order of the singularity at t = 1 increases 
with dimensionality. To approximate r(t), Obreschkow 
et al. make an ansatz based on the general behavior of 
the collapse: 

To = (1 - ?)" (17) 

The exponent a is determined based on the asymptotic 
behavior of the bubble, and is derived in two ways. One 
can match the second derivative of r(t = 0) to r Q (t = 
0), which yields an exponent of t 2 /2 0.418, or one 
can match the behavior of the first derivative at t = 1, 
which yields an exponent of 2/5 = 0.4. Obreschkow et 
al. use the latter, and calculate a series of higher-order 
corrections to find better agreement between r and r a . 

The value of a from the second derivative can be found 
by equating (TT6|) with the second derivative of (fl"7)) . which 
yields the exponent: 

2 4iV r (JV±i) 2 ~ AN [ ' 

For large N the gamma functions reduce to a factor of 
7r. The other method, based on the first derivative at the 
time of collapse, expresses the first derivative in terms of 
the radius: r' a = r\ x l a , This can be equated to (TIT]) 
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FIG. 1: Normalized collapse of a bubble with non- 
dimensionalized radius, pressure, and density with respect to 
time in 3, 6, and 15 dimensions. 

with R a = 1 in the limit R — > 0. This yields the other 
value for the exponent: 



As with the three-dimensional case discussed by 
Obreschkow et al., the first method of calculating a finds 
better agreement with the full Rayleigh collapse equa- 
tion at t=0 while the second method finds better agree- 
ment at t=l. Both methods, however, are within one- 
percent of the numerical solution until the end of the 
collapse. The two values of a are within five percent 
in three dimensions; in the large-iV limit the ratio ap- 

2 

proaches ^ ~ 1.23. The three-dimensional case finds 
the best agreement between the two values, although the 
agreement is exact if the two-dimensional case is allowed. 
Obreschkow et al. chose to focus on the second defini- 
tion of a, and insight from higher dimensions suggests 
that they chose correctly: the dynamics at late times 
in higher dimensions necessitate better accuracy in this 
limit rather than at the onset of collapse. A better un- 
derstanding of the velocity at its zenith is necessary for 
studies of cavitation damage, one of the motivations of 
bubble research. Subsequent to the derivation of their 
approximation, Obreschkow et al. found a series of poly- 
nomial coefficients to improve the precision of their ap- 
proximation, and from the scaling of these coefficients 
found a logarithmic integral that was true to the exact 
result. These coefficients and their scaling can be de- 
rived for arbitrary dimension, but this will be relegated 
to future work. 

The Rayleigh-Plesset Equation 

The Rayleigh-Plesset in three dimensions describes the 
oscillations of a gas-filled spherical bubble. Its derivation 
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is discussed in detail by Leighton (16J, both in terms of 
energy balance methods and the Navier-Stokes equation. 
It is a nonlinear second-order ordinary differential equa- 
tion that takes the following form: 



RR 



3 A. 



-R = 



Pa 



R<> 
~R 



37 



Pc: 



(20) 



The external pressure can include driving terms, such as 
an applied periodic acoustic field, as used in simulations 
of ultrasound contrast agents [17] . The iV-dimensional 
analogue of the Rayleigh-Plesset equation can be derived 
assuming zero mass transfer, demanding the pressure in 
the bubble be a function of its radius, P(R). The work 
performed by the bubble during oscillation is: 

W= / S N R N [P(R) - Pext] dR (21) 
JRo 



and 



By equating the kinetic energy Q and work ([5T 
differentiating, the following relationship emerges: 



Sjyp R 2 R 
~2~W 



N 



f R 

/ S N R N [P(R) - Pext] dR 

JRo 



1], who derives similar expressions in a review of cavi- 
tation science. There is disagreement between the equa- 
tions presented by Prosperetti and those of this paper. 
For example, the iV-dimensional Rayleigh-Plesset equa- 
tion according to Prosperetti merely replaces 3/2, the 
coefficient of the velocity term, with N/2. Obviously, 
experiments cannot be performed in higher dimensions 
to resolve this discrepancy, but the author finds fault in 
Prosperetti's derivation: his equation 3 does not follow 
from his equation 1, as it drops a dimensionally depen- 
dent factor from the second derivative. Other aspects of 
the Rayleigh-Plesset equation in higher dimensions, in- 
cluding the effects of surface tension and viscosity, are 
correctly discussed by Prosperetti. 

There is a case in which the Rayleigh-Plesset equation 
is solvable: that of inertial expansion when the pressure 
term vanishes. Its solution is a power law in time with 
the same exponent as (|19[) . which was also derived by 
examining inertial behavior. 



RR 



NR 2 



N 



2{N - 2) 



-> R(t) = R Q t^ 



(25) 



The Minnaert Frequency 



SnP 
2(N- 2) 



NR 2 R N ~ 1 



2R N R 



S N R N [P(R) - P ext ] 



RR 



NR 2 



P(R) - Pc: 



N 



2(N 



(22) 



Which is the analogue of the Rayleigh-Plesset equation 
in N dimensions. We assume an ideal polytropic gas with 
the following pressure- volume relationship: 



P V = P(P)V 7 P{R) 



P I ^ 
' R 



JV7 



(23) 



Gamma is the polytropic exponent, or ratio of specific 
heats. For an isothermal process, 7 is equal to one. For 
adiabatic systems, it is equal to ^j^, where / is the num- 
ber of degrees of freedom. For a monatomic ideal gas with 
N translational degrees of freedom, f — N and 7 = . 
In the large- N limit, the difference between the adiabatic 
and isothermal cases vanish. Without surface tension or 
viscous damping, and assuming the equilibrium and ex- 
ternal pressures are equal, the full Rayleigh-Plesset ana- 
logue is: 



NR 2 



RR 



N-2 2{N-2) p 



Ro 



if' ' 



(24) 



The Minnaert frequency is the resonance frequency of 
a spherical bubble at low amplitudes. It is derived by 
treating the bubble as a harmonic oscillator and takes 
the following form in three dimensions: 



'3 7 P 

~p~kj 



(26) 



The Minnaert frequency in N dimensions is derived by 
following a similar method as is used for three, again 
assuming small-amplitude oscillations and treating the 
system as a simple harmonic oscillator. The effective 
mass can be calculated by equating the kinetic energy (|9]) 



with the classical definition, 
that of the bubble wall: 

M eS = S N p 



, where the velocity is 



R 



N 



N-2 



(27) 



The effective spring constant can be found by assuming 
R = R + e(t), where e << R , and substituting it into 
the work integral with a first-order Taylor expansion, and 
extracting the spring constant assuming a harmonic po- 
tential: W 



\ke 2 



W 



Jo 



SnR^ 



Ro 



R 



de 



If the equilibrium radius is zero, this becomes a re- 
expression of Rayleigh-like collapse. The dynamics of N- 
dimensional bubbles were also discussed by Prosperetti 
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Evaluating: 



W 



P S N Nj 2 N _ 2 
= e ii n 



Extracting the spring constant we find: 
fceff = S N P Q N 1 R^ 2 



(28) 



(29) 



The Minnaert angular frequency, u>n represents the nat- 
ural frequency of low amplitude oscillations in N dimen- 
sions: 



Voir 



IN(N-2) 1 P 



Mr 



P Rl 



(30) 



The resonance frequency increases with increasing di- 
mension, approaching a linear relationship for high di- 
mension. Like the collapse scenario, the dynamics are 
faster in higher dimensions. One may note that in the 
isothermal case (7 = 1), the Rayleigh collapse time r in 
three dimensions is close to one-quarter of the Minnaert 
period: in dimensionless form, r « 0.914, 2h/Auj « 0.907. 
Simply stated, the time it takes for a perturbed gaseous 
bubble to return to its equilibrium radius is similar to the 
time it takes for a vacuum bubble of the same radius to 
collapse, to within one-percent. This is, however, unique 
to three-dimensions: the ratio of the bubble collapse time 
and perturbation can be found by dividing (| 14[) by (1301) . 
In the large N limit this ratio approaches the square root 
of two. 

Having derived the chief results for bubble collapse and 
oscillation in higher dimensions, the general behavior of 
these systems can be examined. 



DISCUSSION: 



THE ONSET OF NONLINEAR 
BEHAVIOR 



While the total parameter space in which the Rayleigh- 
Plesset equation can be simulated is too vast for a single 
paper, a simple case can be examined to study the onset 
of nonlinear behavior. A dimensionless bubble is stud- 
ied by numerically integrating (|24j) with varying initial 
conditions. 

The Rayleigh-Plesset equation at low amplitudes de- 
scribes a harmonic oscillator. At larger amplutudes, the 
oscillations approach those of a cycloid. In the regime 
R » R a the internal pressure becomes insignificant and 
the dynamics can be described by which is of similar 
form, to the equation of a cycloid. This non-sinusoidal 
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Time (arb) 



FIG. 2: The onset of nonlinear behavior. With larger ampli- 
tudes, the bubble radius spends less time below equilibrium 
and the period begins to shift. 



breathing mode is characterized by a peak wall veloc- 
ity that is much larger near the minimum radius than 
near the maximum, and the oscillations are asymmetric 
about the equilibrium radius, spending much less time 
at small radius (Figure [5]) . In addition to the waveform 
changing from sinusoidal to cycloidal, there is also an up- 
ward shift (in three dimensions) of the oscillation period 
of the bubble, with the shift increasing with amplitude. 
This is analogous to the simple pendulum in the non- 
linear regime, where the period is dependent on the ini- 
tial amplitude. This nonlinear effect has been observed 
experimentally by driving bubble oscillations with ultra- 
sound [18j . 

The onset of nonlinearity was investigated in higher 
dimensions by simulating numerically (|24p with and cal- 
culating the relative radius R(t)/R as a function of 
time. The initial velocity was zero and the initial radius 
was varied from below to above the equilibrium value. 
This defined the independent variable, the initial strain 
R(0)/R o . The nonlinearity can be quantified in several 
ways: the fraction of the oscillation the bubble spends be- 
low its equilibrium radius, the appearance of harmonics 
in a Fast Fourier Transform of the waveform, or the shift 
in the oscillation period or frequency. The period shift 
is primarily is used, as it can be measured from a sim- 
ulation most precisely without dependence on numerical 
parameters and highlights several interesting features. 

The period as a function of initial strain at various di- 
mensions can be seen in Figure [3] In the case where the 
initial strain is greater than one, the oscillation period 
increases with increasing strain in three dimensions. In 
higher dimensions, it decreases and reaches a local mini- 
mum, before increasing monotonically after that minimal 
point. The onset of this behavior as a function of initial 
strain was greater in higher dimensions. For larger ini- 
tial strains, the period increased towards infinity, with a 
greater rate of increase in higher dimensions. The period 
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FIG. 3: The shift in oscillation period as a function of initial 
strain for dimensions 3-10. Dimensions greater than three 
have non-monotonic behavior with a maximal period shift, 
while the nonlinear behavior is stronger at higher dimensions. 

shift approaches a linear relationship with initial radius, 
which is expected when the bubble's oscillations approxi- 
mate Ray leigh- like collapse. The effect on the period was 
symmetric about = 1 for small strains, but asym- 
metry arose due to the fact that the initial size of the 
bubble cannot be less than zero. It is believed that in 
the limit of the initial radius going to zero, the period 
of oscillation doubles. However, computational resources 
were only available to confirm this up to N = 6. This 
limit does not exist for bubbles with initial radii greater 
than R Q . 

The initial strain at which the greatest positive fre- 
quency shift (or negative period shift) occurs was great- 
est in five dimensions (Figure The strain that yielded 
the greatest positive frequency shift occurred approached 
unity with increasing dimension, and the maximal pos- 
itive frequency shift approached the square root of two 
relative to the Minnaert frequency. This ratio, however, 
approaches the square root of two differently than does 
the ratio of Rayleigh collapse time to Minnaert relaxation 
time. 

These simulations also allow the two-dimensional case 
to be examined by considering the limiting behavior of a 
2 + e dimensional bubble. While the period of oscillations 
diverges as e — > 0, the shift in frequency is largely inde- 
pendent when e is small: the frequency shift follows the 
same behavior as the three-dimensional bubble, without 
a local maximum. To quantify the onset of nonlinearity 
as a function of dimensionality, the frequency shift as a 
function of dimensionality for a given initial strain can 
be examined (Figure [5]). The frequency shift at a given 
strain increases monotonically with dimension, with the 
rate of increase dependent on the initial strain. To char- 
acterize the transition from sinusoidal to cycloidal oscil- 
lations, the fraction of the period spent below equilib- 
rium was examined, and showed similar trends to the 
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FIG. 4: The initial strain that leads to a local maximum in 
the oscillation frequency (or minimum in the period, Figure 
fig:bub3), and the maximal frequency at that strain, as a 
function of dimension. It occurs with the greatest strain at 
N — 5, and the relative frequency trends towards the square 
root of two. 

frequency shift with respect to dimensionality: nonlin- 
ear behavior is stronger in higher dimensions. This, as 
well as the dimensional dependence of the period shift, 
highlights the fact that the onset of nonlinearity occurs 
more readily in higher dimension. These simulations 
highlight a few features unique to three dimensions: the 
three-dimensional bubble maintains a sinusoidal breath- 
ing mode at the greatest amplitudes, and the amplitude- 
dependence of frequency is only monotonic. At higher di- 
mensions, there exists a certain large amplitude at which 
the cycloidal oscillations have the same frequency as the 
low-amplitude sinusoidal oscillations, but not in three. 
This behavior is not fully understood and merely serves 
to characterize the nonlinear behavior of bubble oscilla- 
tions in higher dimensions. Because the behavior near 
the minimum size involves wall velocities that approach 
or exceed that of sound in the fluid, a model that takes 
into account compressibility may be necessary to fully 
characterize and understand this behavior. 



CONCLUDING REMARKS 

Several canonical results in bubble dynamics, including 
the Rayleigh collapse time, the Rayleigh-Plesset equa- 
tion, and the Minnaert frequency have been re-derived 
in arbitrary dimension. The bubbles generally display 
faster dynamics in higher dimensions, as the natural pe- 
riod and collapse time both decrease. Particularly, the 
late-time collapse behavior of a bubble is much stronger 
in higher dimensions. The normalized collapse equation 
and its analytic approximation were examined in higher 
dimensions, suggesting a choice of boundary conditions 
upon which to base the approximations. The linear and 
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FIG. 5: Frequency shift as a function of dimension for vari- 
ous initial strains. Nonlinear behavior is stronger at higher 
dimensions. 

nonlinear oscillations of a bubble were simulated numer- 
ically with respect to the initial strain of the bubble's 
radius. It was observed that in higher dimensions the 
onset of nonlinear behavior, characterized by cycloidal 
oscillations with varying frequency, occurs at lesser am- 
plitudes. 

Many factors were not discussed in this paper, includ- 
ing the effect of fluid compressibility and viscosity. Be- 
cause the dynamics in higher dimensions are faster, the 
derivation of a compressible fluid model of ./V-dimensional 
bubble oscillations presents itself as a follow-up. While 
the generalized derivations used in this paper do not ap- 
ply to one- and two-dimensional systems, their quasi-low- 
dimensional real-world analogues have been discussed in 
detail elsewhere [ijj}. Several features were found to be 
particular to three dimensions: the collapse time of an 
empty bubble and the relaxation time of an isothermal 
gas-filled bubble are nearly the same only in three di- 
mensions, and the two dynamical exponents of the ana- 
lytic collapse approximation are closest in three dimen- 
sions. The three-dimensional bubble exhibits nonlinear 
behavior, but this behavior is stronger in higher dimen- 
sions. The shift of frequency with respect to amplitude 
is only monotonic in three dimensions, while bubbles in 
higher dimensions display non-monotonic behavior, hav- 
ing a range of strains that lead to an upward shift in 
frequency. 

Lacking from this paper are true attempts to use 
the higher-dimensional information to better understand 
three-dimensional bubble behavior or solve the relevant 
equations. Solution schemes to the Rayleigh-Plesset 
equation based on perturbation analysis [20( can be 
examined in the light of the .ZV-dimensional Rayleigh- 
Plesset equation to improve upon those methods. These 
derivations may be coupled to higher dimensional models 
of chemical or plasma physics to better understand sono- 
luminescence. It is hoped that in the future more insight 



into three-dimensional bubble dynamics can be gained 
by examining their behavior in higher dimensions. 
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